Can Electrophysiological Parameters Substitute for Growth, and Photosynthetic Parameters to Characterize the Response of Mulberry and Paper Mulberry to Drought?

Drought is a key factor restricting plant survival, growth and development. The physiological parameters of plants are commonly used to determine the water status, in order to irrigate appropriately and save water. In this study, mulberry (Morus alba L.) and paper mulberry (Broussonetia papyrifera (L.) Vent.) seedlings were used as experimental materials, and four soil moisture treatments were set up for both plant species: 70–75% (CK: the control group, referred to as T0), 55–60% (T1: mild drought), 40–45% (T2: moderate drought), and 25–30% (T3: severe drought). The growth parameter of the plants was measured every two days from the onset of the treatment, the photosynthetic and electrophysiological parameters of the plants were measured every other week for a total of five times. The physiological responses and electrophysiological traits of leaves under different treatment levels were analyzed. The results showed that the photosynthetic and electrophysiological parameters could characterize the response of mulberry growth and development to soil water, and the growth and electrophysiological parameters could characterize the response of paper mulberry growth and development to soil water. Mild drought had no significant effects on the growth and development of mulberry and paper mulberry.


Introduction
Drought inhibits plant growth, thereby affecting agricultural production, and its global severity is great concern [1]. With the continuous intensification of the greenhouse effect, the frequency of drought is increasing [2]. Moreover, drought causes morphological and physiological changes, and even causes death [3]. Many parts of China, a large agricultural country, are currently suffering from drought, which not only restricts the development of agriculture and forestry, but also has adverse effects on the ecological environment [4]. Therefore, research on the drought resistance of plants has become a top priority.
Plants can respond differently to drought through complex regulatory mechanisms, which can be summarized as follows: (1) They can enhance water harvesting to escape drought, such as by developing roots, closing stomata, or lowering leaf temperature [5].
(2) They can respond to drought by reducing transpiration, through processes such as premature leaf shedding or accelerated leaf senescence [6]. (3) They can move prematurely from the vegetative to reproductive stage to speed up maturation and seed production before drought induces mortality [7]. (4) They can enhance their drought resistance by improving osmotic regulation and leaf tissue elasticity [8]. (5) The capacity for antioxidant metabolism can be increased to maintain normal growth under severe drought [9] and outside changes, leading to changes in C p , R, and Z in the leaves, and the water inside the cells also changes. The passive electrical properties of plant leaves (C p , R, and Z) can reflect plant water metabolism [34]. Many scholars have explored the response of plant electrophysiology through potential signals, resulting in vague biological significance. Studies of plant electrophysiology were most common until the late 1950s [35][36][37]. Electrical, hydrodynamic, and chemical signal transduction parameters also appear to be regulated when mimosa leaves are subjected to injury stress [38]. The transmission of electrical signals has an effect on photosynthesis and elicits a specific response in leaf photosynthesis [39][40][41]. Studies have shown a significant correlation between the impedance value of leaves and the relative water content, which can be used to accurately evaluate the water content of plants [42].
In recent years, plant electrical signals have been widely used in the field of plant stress resistance physiology in China. Li and Mao used impedance and capacitance to monitor the moisture content of tomato leaves in real time [43]. This research group has also successfully used electrical indicators to detect water status. Zhang et al. defined leaf tension and developed a model of the relationships between leaf tension, tissue water potential and physiological capacitance, which could be successfully applied to rapidly acquire water requirement information in Brassica napus during drought [44]. Xing et al. used the relationship between leaf tension and dry weight biomass to rapidly predict the rehydration time points of leaves under different drought treatments through online monitoring of electrophysiological parameters [45]. We also established a coupling model of leaf clamping force and physiological impedance and found that these parameters play an important role in photosynthesis and water use efficiency. In addition, the use of physiological impedance to monitor water status in plants is nondestructive and can determine the water demand of plants in a timely manner [46]. Since plant cell volume is closely related to physiological capacitance, impedance, resistance and capacitive resistance, electrophysiological characteristics are increasingly used in the diagnosis of plant water status [29,44]. Recently, many cutting-edge technologies have been applied to plant electrophysiology, and an increasing number of studies have been conducted on gene expression related to electrical signals [47].
Mulberry (Morus alba L., M. alba) and paper mulberry (Broussonetia papyrifera (L.) Vent., B. papyrifera), perennial tree species belonging to Moraceae, are characterized by a higher growth rate and greater adaptability to adverse environments than other species in this family [48]. They are widely valuable natural resources. The species M. alba., which is native to China and is now cultivated throughout the world [49], is also an economically important perennial tree, as it serves as the sole food source of the domesticated silkworm. In addition, it adds value through the production of edible fruit, timber and several pharmaceutically important chemicals [50]. B. papyrifera, which is a fast-growing tree mainly distributed in Asian and Pacific countries, can be utilized in many applications. It is well known for its bark fibers, which are used to manufacture high-quality paper, cloth, ropes, animal feed, bioenergy feedstock and traditional Chinese medicines [51,52]. Based on previous studies [53], the two plants have great differences in drought resistance and different mechanisms of water use efficiency. This paper can be used as a comparison material to further study the different characteristics of their electrophysiological responses to drought.
The main contents of this paper are as follows: First, we study the responses of the electrophysiological, photosynthetic and growth parameters of the two plant species to different water treatments. Second, we study the relationships between the electrophysiological, photosynthetic and growth parameters. Finally, we study the different mechanisms underlying the response to drought between the two plant species. Based on these studies, this paper mainly aims to explore the feasibility of using electrophysiological parameters to characterize plant water metabolism, and to clarify the response mechanism of two plant species to drought by combining electrophysiological, photosynthetic and growth parameters.

Soil Moisture Content under Different Drought Levels
As illustrated in Table 1, the moisture content of the soil in which M. alba and B. papyrifera were grown varied depending on the stress level. At the level of T0, the soil moisture content of M. alba and B. papyrifera was the maximum, and it decreased gradually in subsequent levels. At the level of T3, the soil moisture content of M. alba and B. papyrifera was the minimum.
, where m f (g) and m d (g) are the fresh weight of soil and the dry weight of soil, respectively. T0 is the control group, T1 is the mild drought, T2 is the moderate drought, T3 is the severe drought.

Equation-Derived Growth Differences of M. alba and B. papyrifera under Different Water Stress Levels
As illustrated in Tables 2 and 3, the 4-parameter logistic equation fit the plant height, leaf length and leaf width data of M. alba and B. papyrifera under the different water stress levels well. As shown in Tables A1 and A2 in the Appendix A, the coefficient of determination (R 2 ) of all fitted equations is greater than 0.99, and p < 0.0001, indicating that the 4-parameter logistic equation can characterize the relationships between plant height, leaf length and leaf width over time. Based on the GR 50 and DT log values estimated by the 4-parameter logistic equation, the plant height, leaf length, leaf width, and leaf area growth of M. alba and B. papyrifera under drought conditions were calculated (Tables 4 and 5). Note: The means ± SE (n = 3) in the table indicate significant differences in M. alba growth parameters estimated during different water stress phases at p ≤ 0.05, according to one-way ANOVA and Tukey-test. Different lowercase letters in the same column indicate significant differences in the measured factor among the four levels. a is the upper limit of the growth index of the whole growth process; X 0 is the time it takes to reach 50% of the maximum increase in the logarithmic growth period (days); GR 50 is the growth rate in the half of the logarithmic phase; DT log is the duration of the logarithmic phase; and DT s is the time from observation to logarithmic growth. ns is no significant differences. T0 is the control group, T1 is the mild drought, T2 is the moderate drought, T3 is the severe drought. Note: The means ± SE (n = 3) in the table indicate significant differences in B. papyrifera growth parameters estimated during different water stress phases at p ≤ 0.05, according to one-way ANOVA and Tukey-test. Different lowercase letters in the same column indicate significant differences in the measured factor among the four levels. a is the upper limit of the growth index of the whole growth process; X 0 is the time it takes to reach 50% of the maximum increase in the logarithmic growth period (days); GR 50 is the growth rate in the half of the logarithmic phase; DT log is the duration of the logarithmic phase; and DT s is the time from observation to logarithmic growth. ns is no significant differences. T0 is the control group, T1 is the mild drought, T2 is the moderate drought, T3 is the severe drought.  As shown in Table 4, the plant height of B. papyrifera at the T0 level was significantly higher than that at other levels, but the differences were not significant. The leaf length of B. papyrifera at the T0 and T1 levels was significantly greater than that at the T2 and T3 levels, and it was minimal at the T3 level. The leaf width of B. papyrifera at the T0 and T1 levels was significantly higher than that at the T2 and T3 levels. Therefore, the leaf area of B. papyrifera at the T0 and T1 levels was significantly greater than that at the T2 and T3 level, that at the T0 level was the maximal, and that the T3 level was the minimal. Compared with that in T0, the GR 50 × DT log of B. papyrifera in T1 decreased by 7.68%, T2 decreased by 51.57% and T3 decreased by 59.67%. However, as shown in Table 5, the plant height, leaf length, leaf width and leaf area of M. alba were not significantly different among the water stress treatments. The net photosynthetic rate of the leaves of B. papyrifera. T0 is the control group, T1 is the mild drought, T2 is the moderate drought, T3 is the severe drought.

The Transpiration Rate (Tr)
The transpiration rate of M. alba and B. papyrifera showed differences among the water stress levels ( Figure 2). The transpiration rates of M. alba showed obvious differences, with the highest value observed at the T0 level, progressively lower values observed at the T1 and T2 levels, and the lowest value observed at the T3 level, except on 29 June. This pattern was consistent with that observed for the AVG, where the transpiration rate of M. alba was also maximal at T0 and minimal at T3. However, on 1 June, 8 June and 22 June, the transpiration rate of B. papyrifera at the T1 level was significantly higher than that at the other levels. Moreover, the transpiration rate of B. papyrifera at the T3 level was the lowest on 8 June, 15 June and 22 June, but on 29 June, that at the T1 level was the lowest, and that at the T2 level was the highest. The net photosynthetic rate of M. alba at the T0 level was significantly higher than that at other levels at all drought times, that at the T3 level was the lowest, and the AVG showed the same characteristics. On 1 June, 8 June and 15 June, the net photosynthetic rate of B. papyrifera was significantly higher at the T1 level than at the other levels, and that at the T3 level was minimal on 1 June and 8 June. No significant differences were observed on 22 June and 29 June. There was no significant difference in the AVG among the 4-levels.

The Transpiration Rate (Tr)
The transpiration rate of M. alba and B. papyrifera showed differences among the water stress levels ( Figure 2). The transpiration rates of M. alba showed obvious differences, with the highest value observed at the T0 level, progressively lower values observed at the T1 and T2 levels, and the lowest value observed at the T3 level, except on 29 June. This pattern was consistent with that observed for the AVG, where the transpiration rate of M. alba was also maximal at T0 and minimal at T3. However, on 1 June, 8 June and 22 June, the transpiration rate of B. papyrifera at the T1 level was significantly higher than that at the other levels. Moreover, the transpiration rate of B. papyrifera at the T3 level was the lowest on 8 June, 15 June and 22 June, but on 29 June, that at the T1 level was the lowest, and that at the T2 level was the highest. The transpiration rate of the leaves of B. papyrifera. T0 is the control group, T1 is the mild drought, T2 is the moderate drought, T3 is the severe drought.

Stomatal Conductance (gs)
The stomatal conductance of M. alba and B. papyrifera under different water stress levels is presented in Figure 3. The stomatal conductance of M. alba was highest at the T0 level, followed by the T1, T2, and then T3 level, thus showing obvious regularity. The same pattern was observed for the AVG. However, on 29 June, the stomatal conductance of B. papyrifera at the T1 level was higher than that at all other levels, on 8 June, 15 June and 22 June, that at the T3 level was the lowest. On 29 June, the stomatal conductance of B. papyrifera was lowest at the T1 level and highest at the T2 level. The transpiration rate of the leaves of B. papyrifera. T0 is the control group, T1 is the mild drought, T2 is the moderate drought, T3 is the severe drought.

Stomatal Conductance (g s )
The stomatal conductance of M. alba and B. papyrifera under different water stress levels is presented in Figure 3. The stomatal conductance of M. alba was highest at the T0 level, followed by the T1, T2, and then T3 level, thus showing obvious regularity. The same pattern was observed for the AVG. However, on 29 June, the stomatal conductance of B. papyrifera at the T1 level was higher than that at all other levels, on 8 June, 15 June and 22 June, that at the T3 level was the lowest. On 29 June, the stomatal conductance of B. papyrifera was lowest at the T1 level and highest at the T2 level.

The Intercellular CO 2 Concentration (C i )
The intercellular CO 2 concentration of M. alba and B. papyrifera under different water stress levels is shown in Figure 4. The intercellular CO 2 concentration of B. papyrifera showed obvious regularity. On 1 June, that at the T1 level was slightly higher than that at the T0 level, while at other treatment times, the concentration was highest at the T0 level and sequentially lower at the T1-T3 levels, with significant differences observed between the water stress levels. The AVG showed the same characteristics. However, on 1 June, the intercellular CO 2 concentration of B. papyrifera at the T1 level was higher than that at the other levels, and that at the T3 level was the lowest. In terms of the intercellular CO 2 concentration of B. papyrifera, the water stress levels were ordered as T0 > T1 > T2 > T3 on

The Intercellular CO2 Concentration (Ci)
The intercellular CO2 concentration of M. alba and B. papyrifera under differe stress levels is shown in Figure 4. The intercellular CO2 concentration of B. p showed obvious regularity. On 1 June, that at the T1 level was slightly higher tha the T0 level, while at other treatment times, the concentration was highest at the and sequentially lower at the T1-T3 levels, with significant differences observed the water stress levels. The AVG showed the same characteristics. However, on the intercellular CO2 concentration of B. papyrifera at the T1 level was higher tha the other levels, and that at the T3 level was the lowest. In terms of the intercellu concentration of B. papyrifera, the water stress levels were ordered as T0 > T1 > T2 8 June and 22 June, but on 15 June and 29 June, the value at T2 was the highest, at T3 was the lowest. Different letters indicate significant differences at T0-T3 levels in each measurement. Bars with different letters are significantly different at p ≤ 0.05 (Tukey). AVG is the mean value of five measured data. M. alba-g s : The stomatal conductance of the leaves of M. alba. B. papyrifera-g s : The stomatal conductance of the leaves of B. papyrifera. T0 is the control group, T1 is the mild drought, T2 is the moderate drought, T3 is the severe drought.

Water Use Efficiency (WUE)
The water use efficiency of M. alba and B. papyrifera under different water stress levels is shown in Figure 5. The water use efficiency of M. alba showed the same regularity on 8 June, 15 June and 29 June measurement, where it was highest at T0 and decreased sequentially from T1-T3. These differences were significant, and the AVG showed the same characteristics. On 22 June, the concentration at T1 was slightly higher than that at T2, but the difference was not significant. On 1 June, the water use efficiency of M. alba was highest at the T3 level and lowest at the T1 level, but there were no significant differences between the T0, T1 and T2 levels. The AVG showed the same characteristics. At each treatment time, the water use efficiency of B. papyrifera was highest at the T3 level. On 8 June and 22 June, the water use efficiency of B. papyrifera was lowest at T0, and highest at T3. On

Water Use Efficiency (WUE)
The water use efficiency of M. alba and B. papyrifera under different water stress levels is shown in Figure 5. The water use efficiency of M. alba showed the same regularity on 8 June, 15 June and 29 June measurement, where it was highest at T0 and decreased sequentially from T1-T3. These differences were significant, and the AVG showed the same characteristics. On 22 June, the concentration at T1 was slightly higher than that at T2, but the difference was not significant. On 1 June, the water use efficiency of M. alba was highest at the T3 level and lowest at the T1 level, but there were no significant differences between the T0, T1 and T2 levels. The AVG showed the same characteristics. At each treatment time, the water use efficiency of B. papyrifera was highest at the T3 level. On 8 June and 22 June, the water use efficiency of B. papyrifera was lowest at T0, and highest at T3. On 1 June, it was lowest at the T3 level, and on 15 June and 29 June, it was lowest at T2 and highest at T3.  As illustrated in Figure 6, the physiological capacitance of M. alba and B. papyrifera varied with water stress level. The physiological capacitance of M. alba at different drought times showed a decreasing trend across the four water stress levels. The AVG showed the same characteristics. On 8 June and 15 June, although the physiological capacitance of M. alba at the T3 level was higher than that at the T2 level, but the AVG showed the decreasing trend. The physiological capacitance of B. papyrifera at different drought times also showed a decreasing trend across the four water stress levels. On 8 June, 22 June and 29 June, although the physiological capacitance of B. papyrifera at the T1 level was higher than that at the T0 level, but the AVG showed the decreasing trend.

Physiological Resistance (R)
As illustrated in Figure 7, the physiological resistance of M. alba and B. papyrifera varied with water stress level. The physiological resistance of M. alba at different drought times showed increasing trend across the four water stress levels. On 8 June, 15 June and 22 June, although the physiological resistance of M. alba at the T2 level was higher than that at the T3 level, but the AVG showed the increasing trend. The physiological resistance of B. papyrifera at different drought times showed an increasing trend across the four water stress levels. On 8 June and 22 June, although the physiological resistance of B. papyrifera at the T0 level was higher than that the T1 level, but the AVG showed the increasing trend. The water use efficiency of the leaves of B. papyrifera. T0 is the control group, T1 is the mild drought, T2 is the moderate drought, T3 is the severe drought.

Electrophysiological Parameters of M. alba and B. papyrifera under Different Water Stress Levels
As illustrated in Figure 6, the physiological capacitance of M. alba and B. papyrifera varied with water stress level. The physiological capacitance of M. alba at different drought times showed a decreasing trend across the four water stress levels. The AVG showed the same characteristics. On 8 June and 15 June, although the physiological capacitance of M. alba at the T3 level was higher than that at the T2 level, but the AVG showed the decreasing trend. The physiological capacitance of B. papyrifera at different drought times also showed a decreasing trend across the four water stress levels. On 8 June, 22 June and 29 June, although the physiological capacitance of B. papyrifera at the T1 level was higher than that at the T0 level, but the AVG showed the decreasing trend. The water use efficiency of the leaves of B. papyrifera. T0 is the control group, T1 is the mild drought, T2 is the moderate drought, T3 is the severe drought.

Physiological Impedance (Z)
As illustrated in Figure 8, the physiological impedance of M. alba and B. papyrifera varied with water stress level. The physiological impedance of M. alba at different drought times showed an increasing trend across the four water stress levels. On 8 June and 15 June, although the physiological impedance of M. alba at the T2 level was higher than that at the T3 level, but the AVG showed the increasing trend, and the difference of T2 and T3 levels were not significant. However, on 22 June, the physiological impedance of M. alba at the T2 level was obviously higher than that at the other levels, but the AVG showed the increasing trend. In the same way, the physiological impedance of B. papyrifera at different drought times showed an increasing trend across the four water stress levels, with that observed at the T3 level significantly higher than that observed at the other levels. On 8 June, 22 June and 29 June, although the physiological impedance of B. papyrifera at the T0 level was higher than that at the T1 level, the difference was not significant.

Physiological Resistance (R)
As illustrated in Figure 7, the physiological resistance of M. alba and B. papyrifera varied with water stress level. The physiological resistance of M. alba at different drought times showed increasing trend across the four water stress levels. On 8 June, 15 June and 22 June, although the physiological resistance of M. alba at the T2 level was higher than that at the T3 level, but the AVG showed the increasing trend. The physiological resistance of B. papyrifera at different drought times showed an increasing trend across the four water stress levels. On 8 June and 22 June, although the physiological resistance of B. papyrifera at the T0 level was higher than that the T1 level, but the AVG showed the increasing trend.

Physiological Impedance (Z)
As illustrated in Figure 8, the physiological impedance of M. alba and B. papyrifera varied with water stress level. The physiological impedance of M. alba at different drought times showed an increasing trend across the four water stress levels. On 8 June and 15 June, although the physiological impedance of M. alba at the T2 level was higher than that at the T3 level, but the AVG showed the increasing trend, and the difference of T2 and T3 levels were not significant. However, on 22 June, the physiological impedance of M. alba at the T2 level was obviously higher than that at the other levels, but the AVG showed the increasing trend. In the same way, the physiological impedance of B. papyrifera at different drought times showed an increasing trend across the four water stress levels, with that observed at the T3 level significantly higher than that observed at the other levels. On 8 June, 22 June and 29 June, although the physiological impedance of B. papyrifera at the T0 level was higher than that at the T1 level, the difference was not significant.

Discussion
Almost all life activities in plants involve water. Plant physiological information can be represented by the leaf physical parameters which can be directly measured and ensure the accuracy and improve the convenience of drought resistance research of plants [54]. The determination of electrophysiological indexes provides a more convenient method for the study of water status in cells. Electrical signals of plants are considered to be the fastest response of plants to environmental changes and are closely related to the life activities of plants such as material and energy metabolism, development, stress resistance and signal transduction [27,30]. When plant cells lose water and contract, resulting in smaller cell size, the cell expansion pressure must change. Capacitance is related to the degree of expansion and contraction of plant cells [55]. At the same time, the change of water content will inevitably lead to the change of membrane permeability and ion concentration inside and outside the cell, which is closely related to the resistance and impedance [56]. When plants are subjected to different degrees of environmental stress, the water status, ion concentration and membrane permeability of cells will change immediately, and the electrical signals of plants will also change. Therefore, the electrophysiological indexes have the theoretical basis to reflect the water status of plant leaves [57,58]. In this paper, the response of M. alba to drought can be characterized by photosynthetic parameters rather than growth parameters. Drought can cause a decrease in plant cell turgor pressure, which leads to a decrease in components of the plant growth rate, such as plant height, stem thickness and leaf area [59]. Drought has a strong impact on plant growth and development and causes death in severe cases. In addition, the response of B. papyrifera to drought can be characterized by growth parameters rather than photosynthetic parameters. Drought will affect plant photosynthesis, leaf photosynthesis is an important physiological activity of plants, and the accumulated products of photosynthesis play an important role in plant growth and development [60]. However, electrophysiological parameters can characterize the response of M. alba and B. papyrifera to drought. The degree of drought in plants is related to the electrophysiological characteristics of leaves. Therefore, electrophysiological parameters can be used to quickly and nondestructively determine the water status of leaves.
To characterize the effect of drought, the growth curves were analyzed initially using the logistic model, the adequacy of which is discussed. Under different water conditions, the growth models of different plant species are also different [61]. The sigmoidal model measures the growth of plants over time and thus can be used to analyze the temporal growth trend of plants [62]. In this study, a 4-parameter logistic model was used to fit the growth parameters of M. alba and B. papyrifera, such as X 0 , a, GR 50 , DT log , DT s , and R 2 , under different drought levels, and the physiological parameters plant height, leaf length, leaf width and leaf area were analyzed. The results obtained by the 4-parameter logistic model in this study revealed obvious patterns in the growth parameters of B. papyrifera under different drought levels. Changes in the plant height, leaf length, leaf width and leaf area parameters of B. papyrifera showed obvious regularity with the decrease in soil moisture. These results reveal specific life phenomena in plants; that is, B. papyrifera grows fast and has a long-life cycle, and its growth is affected by drought at any treatment time. In contrast, the life cycle of M. alba is shorter, its photosynthesis is lower, and its overall growth is slower. Moreover, the logarithmic growth period of M. alba passed quickly, so drought had no effect on its leaf growth.
Drought resistance in plants can be determined through their photosynthetic and physiological characteristics. The main driving force for the growth of plants is photosynthesis, which supplies the energy and carbon required for the biosynthesis of organic compounds necessary for development [63]. Increasing photosynthetic energy use efficiency and enhanced photosynthetic capacity may be the most successful mechanisms for alien species invasion and for adaptability to adverse environments [64]. In response to water deficiency, stomatal closure occurs, which reduces the net photosynthetic rate and transpiration rate and decreases the intercellular CO 2 concentration [65]. In this study, the leaf area of B. papyrifera is larger than M. alba, and the photosynthesis is stronger than M. alba, the photosynthetic parameters of M. alba responded to soil moisture at different time under different drought treatments, different from M. alba, the photosynthetic parameters of B. papyrifera only responded to soil water in the early stage, but did not respond to soil water in the later stage, indicating that B. papyrifera has strong adaptability. Although the area of M. alba leaves was small, the photosynthetic parameters exhibited obvious differences under different drought levels, indicating that M. alba has a more obvious response than B. papyrifera to water. B. papyrifera continues to grow fast and can be regarded as a "pool". It needs to consume energy for growth. The carbohydrates produced via photosynthesis cannot be consumed through growth, thus inhibiting photosynthesis. Therefore, photosynthetic parameters cannot characterize the response of B. papyrifera to drought. M. alba leaves are small and grow slowly, so the pool is small, which affects photosynthesis. Photosynthesis does not affect the growth of individual leaves. The substance produced via photosynthesis are used to promote the growth of new leaves. M. alba leaves are abundant. Therefore, photosynthetic characters can characterize the response of M. alba to drought.
Plants differ in their physiological capacitance response to water stress. Burdon-Sanderson noted the electrical phenomena that accompany irritation of the leaf of Dionaea muscipula [66,67], prompting studies of plant electrical signals. Previous studies showed that direct changes in plant electrophysiological information such as Z, R and C p could directly reflect changes in plant water [68,69]. Jamaludin et al. study showed that there was a significant correlation between the impedance value of leaves and relative water content, and the water content of plants could be accurately evaluated by measuring the impedance [42]. Yang et al. studied the relationship between plants and air heat in the greenhouse and found that the growth status of plants could be determined by measuring the changes of plant leaf temperature, transpiration and surface electrophysiology of plants [70]. By studying wheat, Briggs found that there was a logarithmic relationship between the water content of wheat and its physiological resistance, and the water content of wheat could be calculated indirectly by measuring the physiological resistance [71]. Ksenzhek et al. studied the physiological resistance of the main vein of maize in different directions, and the results showed that the physiological resistance of the leaf gradually decreased with the wilting of the leaf [72]. In addition, the water status data of plant leaves obtained by photosynthesis-transpiration and water potential did not directly reflect the intracellular water status [73,74]. Through the study of scholars, it is found that the monitoring and analysis of the changes of electrical signals and physiological and biochemical aspects of plants under drought stress, and the test of the response of electrophysiological characteristics of plants to drought stress, so as to achieve a comprehensive, timely and accurate characterization of the drought resistance of crops under drought stress. In this study, the leaf electrophysiological parameters (C p , R and Z) of M. alba and B. papyrifera accurately revealed the diversity of intracellular water metabolism in plant leaves and could potentially be applied in the acquisition of plant water information. The results of this study revealed obvious characteristics of the electrophysiological parameters of M. alba and B. papyrifera With a gradual increase in plant water shortage severity, the capacitance value of M. alba and B. papyrifera decreased gradually and the water content of leaves decreased gradually, and with the increase in resistance and impedance, stomata closed, transpiration decreased, leaf water use efficiency increased, and soil water decreased. Therefore, the electrophysiological index can respond to soil water shortages.
In this study, correlation analysis between photosynthetic and electrophysiological parameters of M. alba and B. papyrifera was performed for the first time (Tables 6 and 7). Plant physiological resistance and intercellular CO 2 concentration are strongly negatively correlated in the two species. The physiological resistance of M. alba and B. papyrifera is also strongly negatively correlated with stomatal conductance; the greater the physiological resistance is, the smaller the intercellular CO 2 concentration. This explains why physiological resistance determines stomatal resistance and physiological resistance leads to greater stomatal resistance and smaller intercellular CO 2 concentrations. The osmotic potential of leaf cells will change immediately upon exposure to drought, so the electrolyte concentration will also change, resulting in changes in leaf capacitance, leaf resistance and leaf impedance. The unique relationships between WUE and electrophysiological parameters indicate the water deficit of plants. The physiological capacitance of M. alba and B. papyrifera is negatively correlated with WUE; the higher the physiological capacitance is, the lower the WUE, indicating that as the physiological capacitance increases, the stomata open, stomatal conductance and the transpiration rate increase, the leaf WUE decreases, and the demand for soil moisture is weakened. At the same time, WUE reflects the unique relationship between the physiological impedance and photosynthetic index of the two plants, because the physiological impedance of plants functions to maintain the internal stability of plant cells [42]. The higher the physiological impedance is, the more stomatal closure will reduce transpiration, while WUE will increase. Therefore, the greater the physiological impedance is, the higher the WUE is and the lower the soil water is, and the electrophysiological index can reflect the soil water demand. This confirms the results of the experiment, which showed that the electrophysiological parameters of M. alba reflected the water requirement. Furthermore, because the electrophysiological parameters were strongly correlated with the photosynthetic parameters, the photosynthetic parameters could also reflect the water demand of M. alba. The water requirements of B. papyrifera can also be inferred by the electrophysiological indicators. However, the electrophysiological indexes were not related to the photosynthetic indexes, so the photosynthetic indexes cannot reflect the water demand of plants. The water demand of B. papyrifera can be reflected only by the electrophysiological indexes and growth indexes.
Concluded from the above, plant electrophysiological parameters have more advantages in characterizing plant drought resistance. Therefore, electrophysiological parameters can substitute for growth, and photosynthetic parameters to characterize the response of mulberry and paper mulberry to drought.

Plant Growth and Treatment
The experiment was carried out in the greenhouse of the Key Laboratory of Modern Agricultural Equipment and Technology of the Ministry of Education, School of Agricultural Engineering, Jiangsu University (N 32 • 11 and E 119 • 27 ). Intact seeds of Morus alba L. and Broussonetia papyrifera (L.) Vent. were sown and germinated in pots containing equal amounts of peat and perlite on March 8, and irrigated with moderate amounts of water. The pots were kept in trays with drainage holes. When the seedlings reached the 3-5-leaf stage, those with uniform growth were transplanted into pots and cultured. After 2 months of seedling acclimation, water control by the weighing method was started in the evening on May 30. Four treatments were set up for both plant species to control the soil moisture content at 70-75% (CK: the control group, referred to as T0), 55-60% (T1: mild drought), 40-45% (T2: moderate drought), and 25-30% (T3: severe drought). Each drought level consisted of three replicates. Water control was carried out for each treatment every day. The measuring time was once every other week from the beginning of the treatment, and the samples were taken in five times, that is, measurements began on 1 June, after 7d of drought stress for 8 June, 14 d of drought stress for 15 June, 21 d of drought stress for 22 June, 28 d of drought stress for 29 June. When the growth of the plants was basically consistent, the plant growth parameters were measured every two days, the photosynthetic and electrophysiological parameters of the plants were measured every other week for a total of 5 times, taken 3 leaves to calculate the mean value.

Determination of Soil Water Content
The drying method was used to measure the soil relative water content. Soil samples were collected with a small shovel into an aluminum box, and the soil samples were weighed with an analytical balance (0.0001 g precision), to obtain the fresh weight of the soil (m f ). The samples were dried in a 105 • C oven for 6-8 h to a constant weight. Then, the dry weight of each soil sample (m d ), expressed by the formula ξ(H 2 O) = (m f − m d )/m d × 100%, was determined.

Determination of Growth Parameters
Growth parameters were measured after treatment application, and three data points were recorded each time. The average value was taken as the growth parameter value of the leaves. The measurements taken for growth analysis were plant height, leaf length and leaf width [69]. S-shaped curves are mainly used to describe the natural growth process of animals and plants and are also known as growth curves. The logistic equation with four parameters best described the growth of the plants [48]. The 4-parameter logistic equation is: where Y 0 is the initial logarithmic growth period; a is the upper limit of the growth index of the whole growth process; X 0 is the time it takes to reach 50% of the maximum increase in the logarithmic growth period (days); the growth rate in the half of the logarithmic phase (denoted by GR 50 ) is GR 50 = −ab/4X 0 ; the duration of the logarithmic phase (denoted by DT log ) is DT log = −4X 0 /b; and the time from observation to logarithmic growth (denoted by DT s ) is DT s = X 0 + 2X 0 /b.

Determination of Leaf Photosynthetic Parameters
The net photosynthetic rate (P N , µmol(CO 2 )m −2 s −1 ), stomatal conductance (g s , mol·m −2 s −1 ), transpiration rate (Tr, mmol·m −2 s −1 ), and intercellular CO 2 concentration (C i , µmol·mol −1 ) were measured with an LI-6400 portable gas exchange measurement system (LI-COR, Lincoln, NE, USA) equipped with a blue/red light source [75,76]. The measurement conditions were as follows: light intensity of 800 µmol·s −1 , CO 2 concentration of 400 µmol·mol −1 , environment temperature of 25 ± 2 • C, the leaf temperature was 30 ± 2 • C, relative humidity of air of 55 ± 5%. These parameters were measured from 9 a.m. to 11 a.m. every seven days under drought treatment. Three leaves were taken from each treatment, and five replicates were taken from each leaf. Water use efficiency (WUE) was calculated according to the equation WUE = P N /Tr, where P N is the net photosynthetic rate and Tr is the transpiration rate.

Determination of Leaf Electrophysiological Parameters
Determination of leaf electrophysiological parameters was carried out according to the method described by Xing et al. [58], with some modifications. The instruments ( Figure 9) included a LCR tester (Model 3532-50, Hioki, Nagano, Japan), a PC (ThinkPad 1430, lenovo, Beijing, China), homemade parallel plate capacitor (1), foam plate (2), capacitor circular electrode plate (3), and wire with a diameter of 10 mm (4). The parallel plate capacitor was connected to the LCR tester with wire, and the LCR tester was connected to the PC. The leaves on the new branches of the plant were selected, and the leaves were clamped between two parallel electrode plates, avoiding the main leaf vein, while connecting the LCR tester and keep the position of the two electrode plates still. The capacitance (C p ), resistance (R), and impedance (Z) were measured using the LCR tester, the frequency and voltage used were 3 KHz and 1.5 V, respectively. Three to five evenly placed points were selected on the leaves of M. alba and B. papyrifera, and 10 measurements were taken in each test. That is, each tested leaf yielded 30 to 50 points (the selected parts of each leaf were guaranteed to be the same), and the experimental data were averaged.
by DTlog) is DTlog = −4X0/b; and the time from observation to logarithmic growth (d by DTs) is DTs = X0 + 2X0/b.

Determination of Leaf Photosynthetic Parameters
The net photosynthetic rate (PN, µmol(CO2)m −2 s −1 ), stomatal conducta mol⋅m −2 s −1 ), transpiration rate (Tr, mmol⋅m −2 s −1 ), and intercellular CO2 concentra µmol⋅mol −1 ) were measured with an LI-6400 portable gas exchange measurement (LI-COR, Lincoln, NE, USA) equipped with a blue/red light source [75,7 measurement conditions were as follows: light intensity of 800 µmol·s concentration of 400 µmol·mol −1 , environment temperature of 25 ± 2 °C, t temperature was 30 ± 2 °C, relative humidity of air of 55 ± 5%.These paramete measured from 9 a.m. to 11 a.m. every seven days under drought treatment. Thre were taken from each treatment, and five replicates were taken from each leaf. W efficiency (WUE) was calculated according to the equation WUE = PN/Tr, where P net photosynthetic rate and Tr is the transpiration rate.

Determination of Leaf Electrophysiological Parameters
Determination of leaf electrophysiological parameters was carried out accor the method described by Xing et al. [58], with some modifications. The instr ( Figure 9) included a LCR tester (Model 3532-50, Hioki, Nagano, Japan), a PC (Th 1430, lenovo, Beijing, China), homemade parallel plate capacitor (1), foam p capacitor circular electrode plate (3), and wire with a diameter of 10 mm (4). The plate capacitor was connected to the LCR tester with wire, and the LCR tes connected to the PC. The leaves on the new branches of the plant were selected, leaves were clamped between two parallel electrode plates, avoiding the main le while connecting the LCR tester and keep the position of the two electrode plates s capacitance (Cp), resistance (R), and impedance (Z) were measured using the LCR the frequency and voltage used were 3 KHz and 1.5 V, respectively. Three to five placed points were selected on the leaves of M. alba and B. papyrifera, and 10 measu were taken in each test. That is, each tested leaf yielded 30 to 50 points (the selecte of each leaf were guaranteed to be the same), and the experimental data were ave

Statistical Analysis
Data were analyzed using exploratory data analysis by SPSS software (version 21.0, SPSS Inc., Chicago, IL, USA) and SigmaPlot software (version 10.0, Systat Software Inc., San Jose, CA, USA). The statistical analysis included a one-way analysis of variance (ANOVA), and significant differences between the means were tested using Tukey test at 95% confidence. The data are shown as the means ± SE. Graphs were prepared using Origin Pro. 9.0 (Northampton, MA, USA).

Conclusions
In both M. alba and B. papyrifera, plant electrophysiological information can be used to quickly describe the response of plants to soil water deficit. There were obvious correlations between the electrophysiological and photosynthetic indexes of M. alba, so these indexes can be used together to represent the response of M. alba growth and development to soil water. Except for physiological resistance, the electrophysiological indexes of B. papyrifera did not show obvious correlations with the photosynthetic indexes, but in general, the growth indexes of B. papyrifera were parallel to the electrophysiological indexes. Therefore, the growth and electrophysiological indexes can be combined to characterize the response of B. papyrifera growth and development to soil water. Mild drought had no significant effects on the growth and development of M. alba and B. papyrifera.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The datasets during or analyzed during the current study available from the corresponding author on reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.